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OPTIMUM  INTERCEPTION  OF  A'  liAj^ElSTIC  MISSILE  AT  MODERATE  RANGE 

The  problem  of  intercepting  a  ballistic  missile  optimally 
.  lit  moderate  z^enge  is  dleouased  in  this  report.  The  problems 
of ,  minimum  time  and  minimum  fuel  consumption  are  discussed 
7’  specifically.  A  procedure  is  given  for  determining  the  trajec¬ 
tory  on  a  digital  computer,  and  some  necessary  and  sufficient 
conditions  are  given  for  the  corresponding  optima. 

The  computational  procedure  is  a  method  of  succepslve  , 
approximations.  The  simpler  cases  have  been  checked  out  on  the 
digital  computer  and  further  programming  is  under  way.  In  the 
examples  run,  each  iteration  takes  somewhat  less  than  one  second 
computing  time  on  the  CDC  l604.  From  three  to  six  iterations 
.  are  usually  required  to  determine  a  trajectory,  so  that  the  tra¬ 
jectory  is  determined  in  five  seconds  or  so. 

The  variational  equations  which  are  used  to  determine  the 
trajectory  can  also  be  used  to  determine  the  corrective  thrust 
'  should  the  actual  trajectory  deviate  slightly  from  the  planned 
.  trajectory. 

The  following  simplifying  assumptions  are  made.  The  speed 
c'  of  the  emitted  gases  of  the  rocket  is  assumed  to  be  a  constant 
of  the  rocket.  The  action  takes  place  above  the  sensible  atmos¬ 
phere  so  that  aerodynamic  forces  are  negligible.  '  The  particu¬ 
lar  cases  are  considered  wherein  the  range  of  the  action  is 
small  enough  that  gravity  is  constant  during  the  flight  or  can 
be  approximated  as  a  linear  function  Of  displacement;  the  case 
of  central  body  motion  is  also  discussed  very  briefly.  Much  . 
of  the  paper  is  restricted  to  the  case  of  plane  mot;ion,  where 
'  the  motion  of  the  rocket  and  the  target  lie  in  a.  properly  defined 
plane.  This  Is  for  simplicity  in  programming  and  not  a  restrlc- 
tlon  on  the  method.  Most  of  the  procedure  is  elementary  and 
hence  the  paper  is  somewhat  expository. 


} 


agitational 


2.  Basic  cquutl'j/is 


The  equation  taf  .Tr.otlci'i  e-f  a  rocket  in  ^  g. 
field  and  with  no  outsxae  tcrces  may  be  written  as 

(1)  r  =  g  +  a, 


wViere  r  is  the  position  vector,  g  Is  the  acoeleratloi^  due 
to  gravity,  a  Is  the  acceleration  due  to  thrust,  and  a  dot 
( * )  over  a  variable  Indicates  Its  time  derivative.  We  may 
write 

(2)  a  =  -  c'm/(l  -  m)  e  , 

where  m  Is  the  ratio  of  the  maos  pf  ’fuel  whlph  has  been 
consumed  to  the  Initial  gross  mass  of  the  rocket  an)i  e  Is  a 
unit  vector  In  the  direction  of  the  thrust;. 

A  useful  kinematic  relation  Is  the  following, 

(3)  pa  dt  =  c* 

Jo 

Since  the  fuel  consumed  Is  proportional  to  m  and  ln{l-m) 

Is  a  monotonlc  function  of  m,  conditions  Involving  the  final 
value  of  mass  may  be  rephrased  In  terms  of  the  Integral  of  the  ^ 
acceleration,  subject  to  the  constraints  on  the  size  of  m. 

For  practical  purposes  we  may  take  m  as  being  bounded  above 
‘by  a  constant,  and  bounded  below  by  zero,  apd  It  may  be  chosen 
anywhere  between  or  on  these  bounds  as  long  as  any  fuel  remains. 


m  dt 
1  -  ,m 


=  -0*  Ind'-m)  . 


Integrated  form  of  the  equations  of  motion.  For  the  case 
where  gravity  may  be  considered  constant,  the  equations  of  motion 


may  be  written 


(4) 


X  =  a  cos  p 
y  =  a  sin  p 


for  two-dimensional  motion,  with  the  y  axis  vertical,.  If 
gravity  may  be  satisfactorily  approximated  by  a  linear  function, 
of  the  displacement,  these  become  , 


(5) 


Where  g 


of  coordinates , 


-b^x 


+  a  cos  p 


y  =  +B  y  +  a  sin  p  “  g  , 

Is  the  gravitational  aoceleratj-op  at  the  origin 


g/l^T 


with 


the  distance  from  the 
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origin  of  coordir.'-  uo  t.l':/:;  o  .  ,.rd  B  =  2  b  . 

The  correcppr.di.r-g  i.a.eanient  may  be  expressed  at  any 
time  T  as 


x(T)  =  Xq  +  XqT  +  (T-t)  a  008  p  dt 

y(t)  =70  +  yo'^  "  gT^/2  +  r  (TTt)  a  sin  p  dt 


'  f  T 

•x(T)  =  x^cos  bT  +  *0  Bln  bT  +  J  |  a  sin  b(T-t)  oos  p  dt 
°  F"  ‘  °  ^0 

^y(T)  a  (yQ+-|)cosh  BT  +  (yQ/B)slnh  BT  -  g/B 
B 

fT 

+'(1/B)  a  slnh  B(T-t)  sin  p  dt. 


3.  Extremals. 

The  curves  which  correspond  to  optimum  motion  from  one 
point  to  another  are  called  extremals.  They  may  be  charac¬ 
terized  for  the  case  where  the  velocity  does  not  enter  In  the 
terminal  conditions  (and  If  gravity  Is  linear,  for  T  <rT/2b) . 
as  follows. 

First,  the  thrust  must  be  a  maximum  for  a  properLy  chosen 
Initial  period,  and  zero  thereafter.  This  Is  a  consequence  of 
the  fact  that  the  coefficient  of  a  cos  p  and  a  sin  p  in 
equations  (6)  and  (7)  diminish  with  time, 

Second,  the  direction  of  the  thrust  Is  glyen'by 


tan  p 


tan  p  = 


:i/B)  slnh  B(T-t) 


]./b)  sin  b(T-ti 


accordingly  as  the  equations  of  motion  are  taken  as  (4)  or  (5). 
Equation  (8)  or  (9)  Is  called  the  steering  equation. 

An  elementary  proof  Is  given  In  section  (5)  to 

show  that  these  conditions  are  sufficient  to  describe  an  ex¬ 
tremal,  and  henceforth  only  trajectories  wherein  t^e  thrust 
satisfies  these  conditions  will  be  conslderad? 


V- 


•J  \J  S-  -L 


The  equeiloii- 
(6)  and  ,(7)  x 


'■  t  he  '  1  u-iOii 

'erla''  .y 


r6c>!<-^  ab  any  time 


*  'l 

,/x(T)  =  x»  +  Xf^T  +  I  (T-t)  a  cos  p  dt 

ho)  .  .1  ,  \  Ti 

ly(T)  =  Yq  +  YoT  “  gT^/2  +  a  sin  p  dt  , 

where  tan  p  =  c,  and 


x(T)  =  XqCOS  bT  +  (xQ/b)  sin  bT  + 

-ft,  ' 

+  d/b'^)  ■^(a/R)sln^b(T-t)  dt  , 

■0 

y(T)  =  (yQ+  g/B^)c08h  BT  +(yQ/BTslnhBT  -  g/B' 

rt 

+  (c/B^)  ^(a/R)slnh^B(TTt)  dt  , 

0 

/«  %  2  .  2^  /  m  l.  \  /  _  /T^  \  2.»  _ ^^2T^/m 


where  R  =  /(l/b)^sin^b(T-t)  +  (c/B)^slnh^B(T-t)  ,  and  t^  <  T. 
the  problem  of  determining  an  optimum  trajectory  is  then  reduced 
to  that  of  finding'  c,  t^,  and  T,  since  these  characterize  an 
extremal . 


4.  Differentials  of  extremals. 

Formulas,  are  derived  here  for  the  changes  in  terminal 
values  of  the  coordinates  due  to  small  changes  in  the  parameters 
of  the  extremals. 

Let  us  consider  now  two  extremals,  near  together  in  the  . 
following  sense.  One  extremal  corresponds  .to  a  sef  of  values 
p„tj^,T,  and  the  other  corresponds  to  a  neighboring  set  p+6p, 
ti+6tj^,T+6T,  in  equations  (10).  Then  the  resulting  change  in^ 
the  terminal  values  of  x  and  y  are 


-  i 


I  '1 

+  ,  ;t  • '  ’a  6p, 

•^0 


for  gravity  constant  if  all  Initial  values  are  specified. 

If  the  linear,  approximation  to  the  gravitational  field  is 
used,  the  equations  for  the  variations  are  more  involved.  In 
the  following  equations  let  C,S,Ch,8ti  denote  cos  b(T-t) , 
sin  b(T-t) , cosh  B(T-t) ,slnh  B(T-t) ,  respectively.  Then 

R  =-/(S/b)^+(c  Sh/B)^,  and 


X  o  2  2 

6X  =  [  a(Gos^p  C  +  2  sin  p  cos  p  C  -  c  cos  p  sin  p  (Jii)dt 

Jo 


-  bxQSln  bT  +  XqCos  bT]&T  +  (aS^/Rb^)^ 

X 


(13)A 


-  (l/c) 


fh 


,6Y  =  .  [ 


a  R  sln^'p  cos  p  dt  6c 

,3, 


a(2  cos'^p  sin  p  Oh  +  sln-^p  Ch  -  -sin  p  cos  p  C)dt 

c 


+  (ByQ  +  |)sinh  BT  +  ^QBinh  BT]6T  +  (ac  Sh^/Rp^) 


(1/c^) 


'■  2  2 
a  R  sin  p  cos  p  dt  6c  , 


f  cb  Sh 

where  tan  P  =  ^  ^ 


5.  Fixed  point  In  fixed  time,  with  minimum  fuel  consumption. 

Let  us,  now  consider  the  problem  of  sending  the  rocket  to 
the  prescribed  point  specified  time  ,  T,  with  minimum- 

fuel  consumption.  The  initial  values  are  all  assumed  to  be 
given.  The  problem  Is  that  of  finding  the  values  of  c,tj^  '  so 
that  x(T)  =  Xg,  y(T)-  =  Yg' 

Suppose  we  guess  a  pair  of  values  for  t^  and  c  (or  p) , 

'  and  compute  the  corresponding  trajectory.  The  end  values  will 
/  not  be  correct  generally;  let  the  errors  be 


Ab-JO-\JVll 


'Now  ict 


E  =  X  ->  X, 


?  *  Y  -  Y, 


6X  =  -E 


Y  =  -F 


and  use  the  formulas  developed  for  6X,6Y.  Thei^ 


/  [(T-t)a  cos'  p].  6t,  -  (T~t)  a  sin  p  dtj  6p  -1 

^1  ^0 


V[(T-t)a  sin  p]^  6t,  +  (T~t)  a  cos  p  dt  6p  =  -F 

^1  ^  '0 

for  constant  gravity,  and  for  the  linearized  equations 

(aS^/Rb^).  6t.,  -  ^  aR  sln^p  cos^p  dt  6c  =  -E 

I  °  Jo 

(17)\  r\ 

I(a3fi^/RB^).  6t,  +  aR  sln^p  cos^p  dt  6c  =  -F. 

^1  ^  )0 

'  The  values  of  6p(or  6c)  ,  6tj^  determined  from  this  yield 

new  values  for  p,t, .  Using  these,  we  compute  a  new  trajec- 

2  2' 

tory  and  correct,  repeating  until  +  F-  Is  less  than 

some  preass Ig-ned  value,  a  convergence  criterion.  No  convergence 
problems  were  encountered  and  a  brief  discussion  of  the  results 
Is  given  In  section  9* 


6.  Proof  of  minimum  fuel  consumption.  A  Mayer  reciprocal 
relation. 

The  proof  Is  given  here  that  the  curve  C*  found  in  the' 
previous  section  has  the  lowest  value  for  J^j'^a  dt  for  admis¬ 
sible  paths.  An  admissible  path  Is  one  whereOn  the  bounds  on 
a  are  satisfied  and  the  proper  Initial  and  final  Values  are 
assumed.  A  corollary  result  Is  the  equivalence  between  problems 
of  specified  displacement  with  minimum  fuel  consumption  and  those 


-7- 
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of  specified  fuel  consumption  and  maximum  displacement,  etc. 
knovm  as  Mayer  reciprocal  relations. 

The  proof  will  be  carried  through  for  equations  (5) J  the  . 
proof  for  constant  gravity  Is  given  In  other  reports  already. 
Suppose  we  have  found  c,t^  In  the  previous  section  so  that 
E,P  are  zero.  Denote  by  an  asterisk  (*)  the  quantities  on 
that  path  C*  and  by  capital  letters  the  corresponding  quan¬ 
tities  on  any  other  admissible  path. 

If  In  equations  (7)  the  second  Is  multiplied  by  c  and  the 
two  are  added,  the  resulting  equation  may  be  rewritten,  for  any 
admissible  path 

(18)  (Xg  -  XqCos  bT  -  fOsln  bT)  +  cCXg  -  ,(yQ+  ^)cosh  BT  - 

b  B 

-  ^Oslnh  BT  + 

B  ■  B^ 

.  f  T  ■  , 

=  a[^'cos.  p  +  gSfi  sin  p]dt  . 


Now  everything  on  the  left  side  of  this  equation  Is  completely 
determined,  either  from  the  Initial  conditions  and  final  con¬ 
ditions,  or  from  c  which  has  been  found.  Hence  the  right 
side  Is  also  determined. 

Now  suppose  there  Is  another  admissible  path  for  which  the 
fuel  consumption  does  not  exceed  that  on  C*.  Thpn 


(19) 

'Let 

(20) 


rT  TT  fT 

A  dt ' <  a*dt  =  1  a*dt . 
0  “Jo  '0 

rh 

(a*-A)dt  =  D  . 

0 


Then,  from  (19)  and  (20) 
fT 

'  A  dt  <  D. 


(21) 


J, 


Now  subtract  the  expressions  obtained  from  (18)  for  the  two  paths. 
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r  t. 


(22) 


■  U. 


a*('^S  cos  p*  +  sin  p*)dt 


rT 


AC^5  oob  P  +  sin  P)dt  ■  0 


Ab-JO-UUli 


Now,  from  equation  (9) 

(23)  ■^S  cos  p*  +  gSli  sin  p*  =  R,  ' 

where  ■ 

=  (i)^S^  +(|)W 

as  defined  earlier.  Note  that  R  Is  a  decreasing  function  of 
t  for  bT  <  ir/2,  and  H(T)  =  0.  In  view  of  this  (22)  may  be 
rewritten  ■ 


(a*-A)Rdt 


+ 


Now 


rt. 


^AHCi  -  cos(p*-P)]dt  - 


rT 


ARoos(p*-P)dt 


"'^1  rt 

(a*-A)R  dt  >  R(t,)  ^(a*-A)dt 

Jo  -  Jo 


=  r(t^)D 


r  T 


AR  oos(p*-P)dt  < 


ARdt 


T 

<  R(tj^)  1  A  dt 


<  DR(t^)  , 

so  that  the  sum  of  the  first  and  last  integrals  In  (24)ls  posl<= 
tlve  or  zero.  The  second  Integral  Is  also  clearly  positive  or 
zero.  Inspection  reveals  further  that  bhe  integrals  In  (2'4) 
have  a  sum  which  Is  not  zero  unless  a*  =  A  except  on  a  set  of 
measure  zero.  Hence  the  curve  found  In  the  preceding  section 
leads  to  a  minimum  value  for  the  fuel  consumption  for  that  value 
of  T  and  Is  unique. 


-9- 
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May«r-reolprooal  rslatlont  Actually  the  trajectory  C*  Just 
found  represents  an  optimun:  relation  among  the  set  of  variables 
Xq ,yQ ,Xq ,yQ jX.Y.mg.  For  consider  all  by  one  of  these  as  being 
prescribed  as  on  C*.  For  example,  take  all  but  Xq  as  having 

prescribed  values.  Let  an  admissible  curve  be  one  on  which 

•  • 

^0 ’*0 ’^0 assume  the  values  they  assume  on  C* ,  with  T 
fixed. 

Then  Xq  assumes  Its  minimum  value  on  C*  for  admissible 
curves. 

To  show  this,  let  us  consider  equation  (18)  for  the  two 
curves.  Everything  on  the  left  Is  prescribed  except  the  term 
Involving  Xq.  Hence 

ft,  fT 

(25)  (-xt  +  X«)cos  bT  =  "^a^R  dt  -  ARcos(p*-P)dt 

^  ^  /Q  JO 

Now  the  right  , side-  of  this  equation  Is  positive  unless  A  =  a*, 
by  the  arguments  given  earlier  In  this  section.  Since  0  <  bT 
<  1^/2,  x^  <  Xq  unless  A  =  a*. 

Similarly,  If  c  Is  positive,  yQ  Is  a  minimum  If  all  the 
other  quantities  are  fixed  as  on  C*.  We  see  that  If  all  but 
one  of  the  set  Is  fixed  and  If  c  Is  positive,  then  C*  defines 
a  minimum  for  Xq ,yQ^Q ,yQ ,m2  and  a  maximum  for  X2,y2. 

’  ..  ' 

.  This  Is  the  Mayer-reciprocal  relation  among  the  set  Just 
discussed  or:  C*.  In  some  cases,  T  may  be  Included  In  the  set 
but  generally  It,  may  not  be.  The  Important  concept  In  the  Mayer- 
reciprocal  relation  Is  that  an  extremal  furnishes  not  Just  a 
maximum  or  a  minimum  for  one  variable  but  rather  an  optimum 
relation  among  a  set  of  variables.  Throughout  the  problems 
solved.  In  this  paper,  we  seek  a  curve  such  as  C*  which  satis¬ 
fies  Just  the  specified  conditions. 


10 
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8.  IntercBpt  1 1.'))!  li'r  mininmri'  •  ime  for  Incoming  target. 

The  problem  of  intercepting  a  ballistic  missile  In  minimum 
time  Is  solved  here.  The  equations  used  are  those  for  the  case 
where  gravity  Is  constant. 

'  The  relative  motion  of  the  missile  and  the  target  may  bo 
written!  In  a  suitably  chosen  coordinate  set  as 


+  X^t  - 


(t-T)a  cos  p  dr, 


(25) 


P  =  Yq  Y^t 


+  Z^t 


(t-'r)a  Bln  p  dT , 

Jo 

0. 


The  number  t*  Is  determined  as  follows.  There  Is  a  maximum 
value  for  t^,  and  In  equations  (25),  t*  =  t  If  t  < 

hmax  ''l  =  hmax  ‘  >  hmax'  ‘"e  computing 

problem  Is  to  make  this  decision. 

To  start  the  computing  routine,  we  chose  the  flrs'p  estimate 

of  p  so  that  the  direction  was  about  60°  from  the  Initial  llni?. 

of  sight  from  rocket  to  target.  The  equations  of  motion 

2  2 

then  Integrated  until  E  +F  began  to  Increase.  This  deter¬ 
mined  the  first  estimate  of  T.  Then  the  variational  equations 
(12)  were  used  with  (25) 


, .  r‘i 

6E,  =  vXq  -  \  a  cos  p  dt)6T  + 


(T-t)a  Bln  p  dt  6p 


(26IY 


a  sin  p  dt) 6T 


0^1 


-^0 


(T-t)a  cos  p  dt  6p. 


To  get  new  estimates  of  T,p  set 


and  solve  for  6T,6p.  This  defines  a  new  approximating  trajec¬ 
tory.  .For  this  and  subsequent  iterations,  the  value  of  T  Is 

determined  In  this  fashion,  rather  than  by  seeking  a  minimum 

2  2  2  2 
value  for  E  +F  .  The  computation  Is  terminated  when  Eg+Fg  , 

Is  below  a  prescribed  number. 


II 


The  oct;  .loove  could  i’cpx'':.,ae:;u  s  maximum 

or  a  rainlmura  value  of  T-  there  are  various  ways  to  check, 

•  One  way  Is  to  dimlnls.  by  an  amount  -6t^  and  calculate 

■  the  corresponding  Increment  6T.  If  6T/5t^  <  0,  then  T  Is 

M  l**it  •  ralatlv®  minimum*  Ths  oasc  whsre  ,thd  rookit  OAn 

"dump"  fuel  Is  not  considered.  If  the  target  Is  not  initially 
Incoming,  then  another  method  must  be  used  to  get  the  first 
estimate  of  T. 


8.  Interception  with  minimum  fuel  consumption. 

The  problem  of  Intercepting  a  ballistic  target  with  mini¬ 
mum  fuel  consumption  Is  discussed  here  for  an  Incoming  target. 
As  In  section  7  the  equations  of  motion  used  are  those  for 
constant  gravity. 

For  Interception,  the,  two  equations  must  be  satisfied 


^0  +  XqT  = 


(28) 


llo  +  = 


(T-t)a  cos  p  dt 
JO  , 

(T-t)a  sin  p  dt 


as  In  the  preceding  section.  There  are  three  parameters 
p,tj^,T  and  only  two  equations  so  there  Is  one  degree  of  free¬ 
dom  In  the  choice  of  the  trajectory.  The  third  equation  Is 
obtained  from  the  condition  dt^/dT  =),  0,  as  follows.  Consider 
two  neighboring  Interception  trajectories,  each  satisfying  (28). 
On  the  second  the  values  p,t^,T  are  replaced  by  p+6p,1:^+6tj^, 
T+6T.  Since  both  trajectories  satisfy  equation  (28), 


.'''t 

1^1 


a  cos  p  dt)6T  -  (T-t^)a( t2^)cos  p  6tj^ 


(29) 


V  i-U/l 


I  ^ 

(^0“  a  sin  p  dt)6T  -  (T-t, )a(tT )sln  p  6t,  - 

I'  ^0  rt  111 

(T-t)a  cos  p  dt  6p  =  0 
Jo 
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When  6p  is  3lin]i.'nated.  ,  the  equation  results 

r^i  T^- 

(39)  ,L(ito  -  1  a  cos  p  dticos  p  +  (Y^  - 
.  JQ  ^ 

(T-t^)a(t3^)6t^  =  0  . 


a  sin  p  dt)sin  p]6T 


The  case  where  =  T  is  an  extreme  case  and  will  be  ex¬ 
cluded  from  consideration  here.  For  a  minimum  value  of  t  , 


rt. 

0 


1  .  fh  , 

(31)  (Xq  -  a  cos  p  dt)cos  p  +  (Xq*“  1  a  sin  p  dt)sln  p  =  0, 


In  formal  calculus  of  variations,  this  is. a  transversal  condition 
Since  the  terms  In  parentheses  are  the  components  of  the  relative 
velocity  vector  for  t  >  t^^,  this  condition  may  be  Interpreted 
as  the  condition  that  the  relative  velocity  be  perpendicular 
to  the  direction  of  thrust,  for  t  >  t^.  If  the  target  follows 
a  known  maneuver,  this  condition  may  hold  only  at  time  T. 

Computational'  routine.  Choose  p  so  that  the  direction  of 
thrust  is  perpendicular  to  the  Initial  line' of  sight.  As  before, 
set 


(25) 


E  =  Xq  +  XQt 


F  -  Yq  “ 


0 

rt. 


(t-T)a  cos  p  dr 


(t-T)a  sin  p  dT  » 


and  guess  t^.  T  may  be  either  be  determined  from-  p  ,  or 
because  the  routine  was  set  up  from  the  preceding  problem, 

■we  computed  E^+F^  until  ,lt  began  to  Increase  to  get  the  first, 
estimate  of  T.  , 

From  equation  (30),  a  new  value  for  p  was  obtained,  and 
hence  6p.  Then  from  the  equations 

rt, 


ftE  =  (Xq  - 


a 


cos  p  dt)6T  -  (T-t^)a(t^)cos  p  6t^  ■ 


rt. 


(3  ){ 


(T-t)a  sin  p  dt  6p 


V 
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(32) 


6F  =  (Yo  ‘~  f  sin  p,  dt)6T  -  (T-ti)a(ti) 

sin  p  btx 


-  I  (T~t)a  cos  p  dt  fip  , 


for  the  variations  and  from 

r  6E  =  -E 
(27)  { 

I  6F  =  -F  , 

6ti,6T  are  determined .  The  value  of  T  is  determined 
from  (27)  for  subsequent  iterations,  rather  than  from 
the  condition  that  E®+F®  is  a  minimum .  The  value  for 
p  is  determined  each  time  from  (31)  and  then  6p,  this 
value  being  used  in  (32) to  reduce  it  to  two  equations  In 
two  unknowns  0 


9o  Comments,  other  problemso 

The  above  method  of  solution  applies  generally  to 
problems  wherein  the  velocity  is  not  involved  at  the  ter¬ 
minal  point.  There  is  also  a  restriction  that  T  be  not 
too  large  if  the  linearized  gravitational  acceleration  is 
used.  In  these  cases,  R  is  always  a  decreasing  function 
of  t,  vanishing  for  t  =  T,  and  the  acceleration  is 
then  always  to  be  applied  initially. 

There  is  a  growing  awareness  that  problems  in  cal¬ 
culus  of  variations  require  the  maximization  of  an  in¬ 
tegral,  and  that  this  in  turn  requires  the  maximization 
of  an  integrand.  The  usual  treatment  of  isoperimetric 
problems  (see  Margenau  and  Murphy  [l],  p,  204)  suggests 
this.  In  formal  calculus  of  variations  this  is  done  by 
an  c-6  procedure,  but,  particularly  in  the  case  of 
linear  differential  equations  it  can,  often  be  done  more 
easily  directly.  The  concept  of  mixlmizing  the  integral 
has.  received  considerable  impetus  recently  due  particularly 
to  papers  by  Pontrjagin  [2], [3],  The  author  first  became 
aware  of  it  from  a  paper  by  Emerson  [4],  and  the  paper  of 
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Breakwell[5]  brings  out  this  property  of  extremals.  The  integral 
in  (18)  is  the  functional  of  Pontrjagln,  and  the  determination 
of  CjtijT  effects  the  construction  of  the  functional. 

The  philosophy  of  solution  here  is  that  one  may  obtain  a 
solution  by  simple  considerations  of  this  maximizing  principle,  ‘ 
and  then  verify  that  it  is  the  solution.  It  should  be  observed 
that  the  proof  given  actually  only  deals  with  the  nature  of  the 
extremals;  the  interrelations  between  the  extremals  and  ’the 
manifolds  whereon  the  terminal  point  lies  are  more  complicated, 
and  we  have  only  discussed  conditions  for  a  stationary  value. 


For  those  interested  in  classic  calculus  of  variations, 
the  condition  1  -  cos  (pip)  ^  0  is  the  Weierstrass  condition 
for  the  steering.  The  strict  inequality  is  of  course  never  satls> 
fled  for  all  values  of  P,  since  cos(p'‘-p)  is  periodic,  but  for 
practical  purposes  it  is  always  satisfied  since  values  of  p 
which  differ  by  2ir  are  equivalent. 


Comments  on  computation.  The  integrals  were  expressed  ^s-^ 
differential  equations  and  the  system  was  integrated  by,,a  Runge- 
Kutta  routine.  It  was  selected  for  convenience  in  programming,  ^ 
with  no  attempt  to  attain  speed.  Each  Iteration,  to  obtain  an 
approximating  trajectory,  requires  about  one  second  on  the  CDC 


1604,  The  squared  error  E^+F^  tended  to  diminish  by  a  factor 
of  60-60  each  time.  The  routine  was  assumed  to  have  converged  ‘ 
when  E2+F2  became  less ' than  in  the  examples.  A 

.further  criterion  should  be  cincluded  corresponding  to  equation 
(31)  but  in  the  examples  run,  which  could  be  checked,  it  was  ^ot 
found  to  be  necessary,  probably  due  to  the  fact  that  in  this 


problem,  the  values  for  the  angle  seem  to  "converge  rapidly". 


It  usually  took  three  to  six  seconds  to  determine  a  trajectory. 
Some  trajectories  were  continued  indefinitely  to  obtain  an 
estimate  of  the  accuracy  of  the  routine,.  In  the  examples  run,  ’ 
the  range  was  a  few  hundred  thousand  feet,  the  time  of  flight 
about  forty  seconds  and  the  initial  acceleration  of  the  rocket 
due  to  thrust  about  6  g  ,  The  terminal  squared  error  ,  E*+F* 
reduced  directly  to  about  10“^®(ft‘®)„ 


Vi 


V  V  -C  i 


An  interesting  convergence  problem  was  encountered  in  the 
problem  of  Interception  in  minimum  time„  The  routine  to  determine 
the  trajectory  was  only  for  an  extreme  value  of  T„  It  seemed  ■ 
obvious  that  for  some  values  of  the  first  estimates  for  c  or 
tan  p  ,  A  minimum  value  of  T  would  be  obtained,  and  for 
other  first  estimates,  a  maximum  value  of  To  Consequently 
a  series  of  trajectories  were  run,  with  various  first  estimates 
for  p,  namely  0,Tr/6,Tr/3,Tr/2,2Tr/3,5ir/6,Tro  The  value  of  p 
associated  with  minimum  T  was  about  lo6,  and  with  maximum  T, 

3,2  (radians).  For  the  initial  estimates  p  =  0,Tr/6,ir/3,Tr/2, 

2n-/3,  the  routine  converged  directly  to  the  value  p  =  1,6, 

For  the  initial  value  p  =  v,  the  routine  converged  directly 
to  p  =  3,2,  For  the  initial  estimate  p  =  57r/6,  the  routine 
did  not  appear  to  converge  initially,  but  the  successive 
estimates  of  p  varied  in  a  random  fashion.  Finally,  p  came 
close  to  14,17'  and  converged  to  this  value,  which  is  l,6+4jr. 

The  method  of  solution  is  essentially  a  Newton’s  method,  and 
this  is  a  typical  behavior  when  the  starting  values  are  not 
neai*  enough  to  the  desired  root. 

These  problems  are  among  the  simplest  in  the  calculus  of 
variations  since  the  equations  are  linear  in  the  dependent  varia¬ 
bles.  The  use  of  a  rather  than  the  mass  simplifies  the  dif¬ 
ferential .  equations  further.  It  can  be  shown  in  several  problems 
that  if  the  time  when  the  fuel  is  to  be  used  is  not  specified, 
then  it  is  to  be  used  when  the  magnitude  of  a  vector,,  R  in 
this  case,  exceeds  a  certain  value  (this  value  to  be  determined), 
and  the  fuel  is  to  be  used  at  the  minimum  rate  whenever  R  is 
below  the  value.  It  is  not  known  whether  this  relation  is 
general  or  notj  no  counterexamples  are  known  to  the  authors. 

The  procedures  for  solution,  using  the  adjoint  system,  can 
be  applied  to  non-linear  problems  equally  well,  though  no  simple 
complete  proofs  exist  corresponding  to  the  ones  given  here.  In 
This. case  the  adjoint  system  is  adjoint  to  the  system  of  varia¬ 
tional  equations  for  the  dependent  variables,  rather  than  the 
original  system,  A  report  is  forthcoming  on  this;  see  also  the 
paper  by  Breakwell  [s],  though  he  does  not  treat  the  case  where 
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a  is  linear  in  the  fuel  consumption  rate„  The  steering  equa¬ 
tion  is  probably  first  due  to  Lawden  in  published  reports  (see 

[6]  for  references) o  The  method  of  calculating  differentials 
is  essentially  that  which  Bliss -formulated  in  Ballistics  ( see 

[7] ,  Chapter  V,  for  introductory  theory,  end  p,  126,  for  early 
pertinent  references) „  The  theory  is  from  lectures  given  at 
Boeing  in  July  1960  on  the  applications  of  the  adjoint  system 
in  control o  An  appendix  is  planned,  giving  the  computationa.1 
details  for  this  problemo 
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Erratum 

The  upper  limit,  of  the  last  integral  in  equation  (19)  should  be 
ti,  not  To 
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